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When  turbulent  jets  or  wakes  interact  with  a  free 
surface,  the  turbulent  eddies  are  apparently  damped  in  the 
vertical  direction  and  extended  in  the  horizontal  direction. 
Turbulent  velocities  are  being  measured  for  this  zone  of 
interaction  to  enable  verifying  turbulence  modelling 
techniques.  A  two-dimensional  jet  is  formed  by  pumping 
water  through  a  rectangular  slit  into  a  channel  filled  with 
still  water.  Hot-film  anemometry  techniques  are  used  to 
measure  velocity  fluctuations  in  the  jet  and  in  the 
immediate  neighborhood  outside  the  jet.  Some  effects  of  jet 
flow  rate  and  jet  submergence  have  been  studied.  . 

The  raw  data  are  sampled  time  records  of  the  anemometer 
output  voltages  obtained  at  different  positions  in  the  jet. 
These  time  records  together  with  the  sensor  calibration 
information  are  analyzed  using  digital  signal  processing 
techniques.  The  mean  velocities,  turbulence  intensities, 
and  velocity  correlation  functions  are  (currently  being 
calculated* 

The  results  are  being  used  to  develop  a  two-point 
closure  scheme  and  a  mathematical  model  for "the^  jet  flow 
turbulence.  The  model  contains  small  non-homogeneities  and 
anisotropy  because  of  the  proximity  to  the  free  surface.  By 
analogy,  the  results  are  being  used  to  predict  the 
interaction  of  a  wake  with  a  free  surface. 


Background  and  Objectives 

The  wake  of  a  ship  or  torpedo  gives  evidence  of  its 
passing.  For  detection  and  tracking  purposes,  it  is  of 
interest  to  know  as  much  as  possible  about  the 
characteristics  of  wakes.  It  is  difficult  to  study  wakes 
formed  by  moving  objects,  because  the  detectors  themselves 
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must  be  moving.  A  towing  channel  or  basin  is  required  for 
this ,  and  moving  the  detector  can  introduce  extraneous 
vibrations  which  can  make  data  analysis  difficult.  The  flow 
in  a  wake  can  be  simulated  by  the  flow  in  a  jet.  This  is 
one  unique  feature  of  the  current  research  being  supported 
by  ONR. 

The  similarities  between  turbulent  flow  in  jets  and 
wakes  have  been  described  by  Schlichting(1968) .  The  mean 
velocity  profiles  have  similar  shapes,  but  the  jet  profile 
is  the  negative  of  the  velocity  defect  of  the  wake  profile. 
The  width  of  a  jet  or  wake  is  proportional  to  the  distance 
downstream  from  the  source,  and  the  centerline  velocity  is 
proportional  to  the  square  root  of  the  distance  downstream. 
The  similarity  becomes  clear  if  the  velocity  defect  profile 
for  the  wake  is  compared  to  the  velocity  profile  for  the  jet 
as  shown  by  Abramovich(1963) .  since  the  characteristics  of 
jets  and  wakes  are  similar,  studying  jet  turbulence  should 
give  insight  into  wake  turbulence. 

In  the  current  research,  the  interaction  of  a  two 
dimensional  plane  jet  with  a  free  surface  is  being  studied 
by  measuring  the  distribution  of  fluctuating  velocities  in 
the  flow  from  a  submerged  jet  located  near  the  surface. 

This  is  being  done  by  making  time  records  of  the  output  of 
hot-film  anemometers  placed  at  various  positions  in  the  jet. 
From  these  time  records  and  sensor  calibration  data,  the 
mean  velocity  distribution,  jet  width,  turbulence  intensity, 
and  turbulent  velocity  correlations  are  being  calculated.  A 
turbulence  model  based  on  a  two-point  closure  scheme  is 
being  developed,  and  the  results  of  these  modelling 
calculations  is  being  compared  with  the  experimental  data 
and  with  data  which  has  been  reported  by  others. 

The  objectives  of  the  current  project  are  as  follows. 

(1)  Determine  the  effect  of  jet  flow  rate,  depth  of 
submergence ,  angle  of  impingement,  and  jet 
dimensions  on  ther  turbulence  characteristics 

in  the  zone  where  a  two-dimensional  turbulent  jet 
interacts  with  a  free  surface.  These  turbulence 
characteristics  include— 

-  steady-state  mean  flow  velocity 
distribution 

-  distribution  of  turbulent  intensity 

-  spatial  distribution  of  turbulent  velocity 
correlations 

-  time  decay  characteristics  of  velocity 
fluctuations  and  velocity  correlations 

(2)  Develop  a  two-point  closure  model  for  turbulence 
near  a  free  surface.  This  model  includes  some 
interesting  effects. 

-  Small  nonhomogeneities  occur  because  of  the 


presence  of  the  free  surface 
-  Anisotropy  is  present,  because  there  is  a 
preferred  direction  perpendicular  to  the 
free  surface. 

(3)  Determine  the  time  decay  properties  of  turbulent 
velocity  correlations  as  a  consequence  of  the 
closure  scheme. 

The  literature  pertinent  to  this  project  was  reviewed 
in  the  original  proposal.  Most  of  the  papers  dealing  with 
experimental  measurements  describe  studies  made  using  air 
jets  where  there  was  no  anisotropy  caused  by  the  presence  of 
an  air-water  interface.  Many  of  the  mathematical  modelling 
papers  in  the  literature  deal  with  one-point  closure 
models— the  k-epsilon  model  ,for  example.  These  cannot 
incorporate  fully  the  turbulent  energy  transfer  from  low 
wave  number  eddies  to  high  wave  number  eddies.  The 
contemporary  view  of  steady-state  turbulence  is  that  flow 
structures (coherent  structures)  or  eddies  decay  in  the  flow 
and  eventually  disappear.  In  order  to  sustain  a  steady- 
state  turblent  flow,  new  structures  must  be  formed,  and  this 
requires  energy.  The  energy  is  thought  to  come  from  the 
mean  flow,  but  the  mechanism  for  the  transfer  from  the  mean 
flow  to  the  turbulent  structures  is  not  fully  understood. 
WHATEVER  THE  ENERGY  TRANSFER  MECHANISM  IS,  ITS  EFFECTS 
SHOULD  BE  REFLECTED  IN  ANY  CLOSURE  MODEL  DEVELOPED  IN  ORDER 
FOR  THE  MODEL  TO  BE  UNIVERSALLY  APPLICABLE  TO  ALL  TURBULENT 
FLOWS. 

Two-point  turbulence  closure  modelling  techniques  for 
isotropic  turbulence  have  been  discussed  by  Trevino 
(1982a, 1983)  and  by  Domaradski  and  Mellor(1984) .  The 
Navier-Stokes  equations  governing  momentum  transport  in 
fluids  are  written  in  terms  of  the  two-point  velocity 
correlation  tensors. 


|^Ci;J(r,t)  -^Sikj(r,t)  -~Sjki (-£, t)  +  Pressure  terms 
^  +  Viscous  terms 


This  equation  is  being  solved  by  assuming  self-similar 
structure  in  the  turbulence,  that  is,  by  assuming  that  the 
turbulence  and  the  mean  flow  scale  in  the  same  way 
(Trevino, 1986a) ,  and  by  including  the  "free-surface"  effects 
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Accompl ishment s  Purina  The  First  Year 

Work  on  the  jot  flow  research  project  began  on 
1  March  1985.  The  firat  four  months  were  devoted  to 
equipment  design,  construction,  and  installation.  This  was 
a  time  consuming  process,  since  an  entirely  new  project  was 
being  started. 

A  special  device  for  producing  a  plane  jet  was  designed 
and  fabricated.  There  are  entrance  sections  and  a  narrow 
slit  as  shown  in  Figure  1.  The  parallel  flat  plate  section 
forming  the  slit  is  separate  so  that  slit  sections  having 
different  slit  widths  can  be  attached  to  the  entrance 
sections  enabling  the  effect  of  slit  width  to  be 
investigated.  This  jet  flow  unit  is  installed  in  a  water 
channel  0.45  m  wide,  0.45  m  deep,  and  7.5  m  long.  A 
schematic  flow  diagram  of  the  system  is  shown  in  Figure  2.  A 
special  framework  is  installed  over  the  channel  to  support 
the  jet  flow  unit.  This  framework  keeps  the  unit  rigidly  in 
position  and  incorporates  a  means  for  raising,  lowering,  and 
tilting  the  jet  flow  unit.  Water  is  supplied  to  the  jet 
flow  unit  from  a  pump  supported  above  the  channel  on  a 
separate  framework.  The  piping  incorporates  a  flexible 
tubing  section  to  isolate  the  jet  flow  unit  from  pump 
vibrations.  The  pump  suction  line  is  located  downstream 
from  the  jet  to  minimize  reverse  flow  effects. 

The  fluctuating  velocities  in  the  jet  are  measured 
using  hot-film  sensors  connected  to  a  TSI,Inc.  Model 
IFA-100-4  anemometer  system  which  was  purchased  partly  with 
funds  from  ONR  and  partly  with  funds  from  Michigan  Tech. 

This  four-channel  anemometer  system  enables  obtaining  time 
histories  of  the  fluctuating  velocities  simulatneously  at 
four  different  points  in  the  flow.  This  enables  determining 
six  points  on  the  spatial  correlation  function  curve.  This 
gives  a  definitive  test  of  the  degree  to  which  the 
experimental  data  deviate  from  the  corresponding  curve  which 
can  be  calculated  for  isotropic  turbulence  and  enables 
making  a  good  estimate  of  the  inherent  anisotropy  which 
occurs  because  of  the  presence  of  the  free  surface.  The 
hot-film  probes  are  mounted  on  a  framework  separate  from  the 
other  two  to  eliminate  any  vibration  which  might  be 
transmitted  from  the  pump  or  elsewhere.  The  probe  can  be 
positioned  anywhere  in  the  jet  by  sliding  the  mounting 
device  horizontally  along  a  set  of  rails.  At  any  horizontal 
position,  the  probe  support  can  be  moved  to  any  vertical 
position  and  clamped. 

One  problem  which  occurs  when  operating  hot-film  probes 
in  water  is  fouling.  This  cannot  be  avoided,  so  the  probes 
must  be  calibrated  frequently.  It  is  essential  to  have 
timely  and  accurate  calibrations  in  order  to  be  able  to 
interpret  the  anemometer  output  voltage  data.  It  would  be 
very  time  consuming  to  calibrate  each  probe  frequently 
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during  the  course  of  an  experiment.  Because  of  this,  quite 
a  bit  of  time  during  the  summer  of  1985  was  spent  to  develop 
a  new  type  of  laminar  jet  calibration  system  and  a  new 
calibration  procedure.  This  device  uses  water  from  the  same 
source  as  that  which  is  flowing  through  the  slit  which 
produces  the  turbulent  jet.  The  construction  of  the 
calibration  device  is  shown  schematically  in  Figure  3.  It 
is  positioned  as  shown  in  the  flow  diagram (Figure  2).  The 
method  of  using  this  calibration  device  to  incorporate  probe 
fouling  while  making  turbulent  velocity  measurements  is 
discussed  in  a  paper  by  Hubbard,  Trevino,  and  Hine(l985) 
which  was  presented  at  the  22nd  Annual  Meeting  of  the 
Society  of  Engineering  Science.  This  paper  also  includes  a 
discussion  of  all  the  operating  procedures  which  have  been 
developed  in  connection  with  the  jet  flow  research  project. 

A  copy  of  this  paper  is  attached  as  an  appendix.  The 
development  of  the  probe  calibration  scheme  and  the  method 
for  handling  probe  fouling  is  one  of  the  unique 
accomplishments  for  the  project. 

Besides  completing  the  design,  construction,  and 
installation  of  the  equipment  and  the  development  of 
operating  and  calibration  procedures,  computer  software 
necessary  for  carrying  out  the  data  analysis  has  been 
written.  The  anemometer  output  data  and  the  calibration 
data  are  transferred  to  the  mainframe  computer  into  separate 
files.  The  voltage-time  data  are  converted  to  velocity-time 
data  by  applying  the  appropriate  calibration  function.  The 
new  files  thus  formed  contain  the  instantaneous  velocity 
values.  These  data  can  then  be  analyzed  using  the  time¬ 
averaging  codes  or  the  velocity  correlation  codes. 

Each  set  of  voltage-time  data  consists  of  data  obtained 
at  the  same  time  at  four  different  locations  downstream  from 
the  source  of  the  jet.  Using  these  data,  time-averaged 
velocities  have  been  calculated  to  show  the  validity  of  the 
original  idea  of  using  the  turbulent  jet  as  a  model  system 
to  represent  a  turbulent  wake  and  to  obtain  some  idea  about 
the  persistence  of  the  jet  with  distance  from  the  source. 

The  results  of  some  of  these  calculations  appear  in  Figure  4 
where  the  mean  velocity  profiles  for  the  jet  are  compared 
with  data  obtained  at  the  Naval  Research  Laboratory  by  Swean 
and  Keramidas(1984) .  These  data  show  that  the  two  systems 
are  indeed  similar,  and  the  original  idea  of  using  jets  to 
model  wakes  is  verified. 
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Correlation  calculations  are  presently  being  made  in 
preparation  for  the  detailed  turbulence  closure  modelling. 

To  make  the  correlation  calculations,  it  is  necessary  first 
to  calculate  the  time-averaged (or  mean)  value  for  the  data 
in  each  file  and  subtract  this  value  from  all  the  velocity 
values  to  obtain  the  turbulent  velocity  fluctuations.  The 
correlation  function  values  or  double  velocity  correlations 
are  calculated  at  six  different  values  of  the  separation, £. 

Cij(r)  -  <  Ui'fs)  Uj  •(*+!)  > 

where  u'  is  the  velocity  fluctuation  and  the  <  >  brackets 
indicate  an  ensemble  average  which  is  approximated  here  by 
the  time-average.  The  skewness  or  triple  velocity 
correlations  are  also  calculated  for  24  different 
combinations  of  the  separations  available  from  the  data. 

sijk<*)  *  <  ui'(X)  Uj'(X)  uk*(x+r)  > 

As  part  of  the  work  on  the  contract,  George  Trevino  has 
developed  the  general  theoretical  framework  for  the 
turbulence  closure  scheme  based  on  self-similar  velocity 
correlations.  A  paper  entitled  "Analysis  of  Turbulent  Decay 
Using  Affine  Transformations"  has  been  prepared  and  is 
currently  being  reviewed  for  publication  in  Physics  of 
Fluids.  Another  paper  entitled  "The  Skewness  Function  of 
Turbulence"  has  been  completed  and  submitted  for 
publication  in  Boundary  Laver  Meteorology.  Copies  of  these 
manuscripts  are  attached  in  the  appendix. 

The  idea  of  self-similar  structure  means  that  the 
turbulence  characteristics  retain  similarities  throughout 
the  turbulent  field.  This  does  not  mean  that  these 
characteristics  must  be  indent ical  throughout  the  field.  It 
only  means  that  something  about  the  turbulence 
characteristics  is  preserved  throughout  the  field.  For 
example,  the  ratio  of  the  integral  scale  to  the  Taylor 
microscale  might  be  preserved  as  the  turbulence  decays— as 
the  distance  downstream  from  the  source  of  the  jet 
increases.  This  ratio  of  scales  is  related  to  the  triple 
velocity  correlation  or  skewness  of  the  probability  density 
function  which  characterizes  the  jet  turbulence. 

In  Figure  5  are  shown  some  triple  correlation  data 
which  illustrate  this  idea  of  self-similar  structure.  These 
data  show  that  the  triple  correlation,  S^jk,  depends  on  the 
separation  distances  but  not  on  the  distance  downstream  from 
the  source  of  the  jet,  because  all  the  functions, S(x,e,r) 
have  the  same  shape  without  regard  to  the  position  in  the 
jet.  Since  just  four  sample  records  are  available,  only  a 
few  points  on  the  triple  velocity  correlation  function  can 
be  determined,  but  the  data  available  are  consistent  with 
the  idea  of  self-similar  decay. 
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The  complete  results  obtained  during  the  first  year  for 
the  jet  flow  project  will  be  presented  in  the  thesis 
currently  being  written  by  Eric  J.  Hine  who  was  supported  by 
funds  from  the  contract  while  pursuing  his  graduate  work.  A 
paper  based  on  this  work  is  being  prepared  and  will  be 
presented  at  the  annual  meeting  of  the  American  Institute  of 
Chemical  Engineers  in  November  1986.  A  copy  of  the  abstract 
of  this  paper  and  a  copy  of  the  program  for  this  session  are 
attached  in  the  appendix.  Eric  Hine  has  accepted  employment 
at  the  Knolls  Atomic  Power  Laboratory  and  will  be  working  on 
nuclear  ship  propulsion  systems.  He  will  begin  in  this 
position  as  soon  as  he  finishes  his  thesis  and  his  security 
clearance  is  completed. 


Plans  for  Continuing  the  Project 

During  the  next  year,  it  is  proposed  to  continue 
turbulent  velocity  measurements  to  provide  more  data  for 
comparison  with  turbulence  models.  These  data  will  be 
obtained  for  different  values  of  jet  submergence  and  for 
different  values  of  the  angle  of  impingement  on  the  free 
surface,  so  that  these  effects  can  be  studied.  It  is  also 
proposed  £2  add  pressure  measurements  to  the  experimental 
program.  The  closure  and  turbulence  modelling  calculations 
will  be  continued,  and  the  model  will  be  improved  by 
incorporating  pressure-velocity  correlations  into  the  model. 

Cheng(1985)  and  Trevino (1986b)  point  out  the  importance 
of  pressure-velocity  correlations  in  modelling  the  spatial 
evolution  of  the  transfer  of  turbulent  energy  to  lower  and 
higher  wave  numbers.  For  isotropic  turbulence,  the 
pressure-velocity  correlation  is  zero,  and  many  turbulence 
closure  models  are  based  on  setting  it  equal  to  zero.  For 
non-isotropic  turbulence,  the  pressure-velocity  correlations 
are  not  zero.  It  is  proposed  to  compute  the  pressure- 
velocity  correlations  and  use  the  results  to  introduce  more 
realism  into  the  closure  model  being  developed  as  part  of 
the  jet  flow  project. 

In  the  "spatial  decay"  fspatlal  evolution)  of  a  gt?.fl<3y- 
state.  self-similar,  nonhomogeneous  turbulence  such  as  may 
be  encountered  in  a  submerged  turbulent  jet,  the  convective 
mechanism  is  spatially  invariant;  that  is,  independent  of 
spatial  location  downstream  within  the  jet.  This  is  true 
even  though  the  turbulence  itself  is  fully  nonhomogeneous. 
This  point  is  discussed  by  Trevino  in  some  forthcoming 
papers(see  Trevino, 1986a, b) .  This  spatial  independence 
therefore  requires  that  in  the  "high"  Reynolds  number 
situation,  it  is  the  pressure-velocity  correlation  which 
provides  the  all-important  spatially  varying  energy  transfer 
mechanism.  Spatial  invariance  of  the  convective  tensor, 

S 4 4 1 ( . . . )  -  <U4  (2£— E/2)U4  (H-I/2)u1(x+e/2)>,  can  exist  only  if 
there  is  a  continuous  cascade  of  energy  toward  the  lower 
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frequencies  of  the  spectrum  as  well  as  toward  the  higher 
frequencies.  This  phenomenon  has  already  been  suggested 
theoretically  by  Tsuge(1984)  and  shown  experimentally  by 
Okude(1981) .  The  role  of  P*(x,r)  -  <P(X-r/2)uA (*+£/2) >  in 
the  decay  of  -  <ui  (2£-£/2)uj  (x+i/2)>  accordingly  ought 
not  be  neglected (  as  it  is  typically  neglected  in  the 
analysis  of  nonhomogeneous  decay) ,  and  its  effects  must  be 
explicitly  measured  and  specifically  included  in  the 
equation  of  jet-decay.  Because  of  the  anisotropy  present  in 
the  turbulence,  the  pressure-velocity  term  is  not  zero.  The 
related  steady-state  dynamical  equation  of  turbulence  for 
high  Reynolds  number  is — 

U  *Viz  +  a  *75  +  5  *7a  +  a  *7a  -  -  £  Vp 

where  the  total  flow  is  !Z  a  Q  +  u  end  the  total  pressure  is 
P  ■  P  +  p.  In  the  process  of  stochastic-averaging,  this 
equation  produces  the  Nspatial-decayN  equation— 

<11-72  a’>  +  2  7<u  a'>  +  <a*7a  a'>  ■  -pV<pa'> 

where  a^a(X“iy2)  and  a(X+l/2) ,  etc.  For  the  self¬ 
similar  turbulence  being  considered,  the  term  <a*Va  a’>  is 
independent  of  j£#  but  the  term  (-  1/f  )7<p  a'>  i»  not. 

This  latter  term  is  in  fact  what  introduces  into  the  right- 
hand  side  of  the  equation  the  ^-dependence  necessary  to 
sustain  the  ^-dependence  of  the  terms  <a*7fl  a'>  end 
5*  7<a  u'>.  THIS  PRESSURE  TERM  IS  THE  TERM  FOR  WHICH 
MEASUREMENT  AND  MODELLING  ARE  CRUCIAL  TO  THE  ANALYSIS. 

Until  recently,  pressure  sensors  sensitive  enough  to 
measure  turbulent  pressure  fluctuations  were  not  available, 
so  the  effect  of  the  pressure  terms  could  not  be  included  in 
two-point  turbulence  closure  models  for  nonisotropic 
turbulence.  Jones  (1981)  and  Spencer(1970)  report  the 
development  of  a  bleed-type  pressure  sensor  based  on  a  hot- 
film  anemometer  system.  A  hot-film  velocity  probe  connected 
to  a  constant  temperature  anemometer  circuit  is  enclosed  in 
a  tube  open  to  the  pressure  to  be  measured  at  the  downstream 
end.  The  upstream  end  is  connected  to  a  constant  pressure 
source  of  fluid,  and  there  is  a  constant  flow  of  fluid  from 
the  constant  pressure  supply,  past  the  velocity  sensor,  and 
into  the  flow  where  the  pressure  is  being  measured.  The 
velocity  of  the  fluid  flowing  through  the  tube  is  related  to 
the  pressure  difference,  so  the  downstream  pressure  can  be 
related  directly  to  the  anemometer  output.  It  is  well-known 
that  hot-film  velocity  probes  are  very  sensitive  to  changes 
in  the  cooling  velocity  in  the  neighborhood  of  the  sensor. 
This  means  that  pressure  measurements  based  on  the  velocity 
determined  using  a  hot-film  sensor  are  equally  as  sensitive. 
Pressure  sensors  such  as  described  above  are  available 
commercially  from  TSI,  Inc.  For  measurements  in  water,  the 
bleed  fluid  used  is  alcohol  which  is  miscible  with  the  main 
water  flow. 
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The  research  program  in  progress  is  yielding  a  body  of 
results  which  will  indicate  the  effects  of  jet  flow  rate, 
jet  submergence,  and  jet  orientation  on  the  turbulence 
characteristics  and  decay  in  the  zone  of  interaction  between 
a  plane  jet  and  a  free  surface.  When  the  jet  is  submerged 
well  below  the  surface,  it  is  expected  that  the  time- 
averaged  velocities,  turbulence  intensities,  Reynolds 
stresses,  and  other  turbulence  characteristics  will  be  the 
same  as  similar  data  reported  by  others  for  air  and  water 
jets.  The  results  for  jets  at  shallow  submergence  are  new 
data  which  can  be  compared  with  theoretical  results  for 
anisotropic  closure  schemes.  The  comparison  will  establish 
the  range  of  validity  of  the  anisotropic  decay  laws 
generated  by  the  theory. 

Besides  the  one-point,  steady-state  turbulence 
characteristics,  the  experimental  results  will  include 
information  about  the  decay  of  the  turbulence  parameters  and 
information  about  the  spatial  dependence  of  the  two-point 
velocity  correlation  functions.  The  two-point  closure  model 
to  be  developed  will  include  pressure-velocity  correlations 
as  well  as  spatial  decay  effects.  Because  of  the  presence 
of  the  free  surface,  small  nonhomogeneities  and  anisotropy 
will  also  be  included.  The  two-point  model  is  more 
effective  for  modelling  the  details  of  the  turbulence  far 
from  the  source  of  a  jet  or  wake. 


The  project  is  being  carried  out  by  D.  W.  Hubbard 
(experimental  work)  and  G.  Trevino  (closure  modelling)  with 
the  help  of  graduate  research  assistants  who  are  candidates 
for  advanced  degrees  in  chemical  engineering  or  mechanical 
engineering.  The  research  is  the  basis  for  preparing  theses 
which  will  fulfill  part  of  the  requirements  for  these 
degrees.  The  students  receive  academic  credit  for  the  work 
and  spend  approximately  one-half  time  engaged  in  the 
research.  The  students  currently  working  on  the  project  are 
Eric  J.  Hine  and  Lisa  Heydenburg  who  are  candidates  for 
graduate  degrees  in  chemical  engineering. 
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EXPERIMENTAL  TECHNIQUES  FOR  MEASURING  TURBULENT  VELOCITY 
CORRELATIONS  IN  A  JET  FLOW  NEAR  AN  AIR-WATER  INTERFACE 


D.  W.  Hubbard,  G.  Trevino,  and  E.  J.  Hine 

Michigan  Technological  University 
Houghton,  Michigan  49931  USA 


ABSTRACT 


When  a  turbulent  water  jet  interacts  with  an  air- 
water  interface,  the  turbulent  eddies  are  apparently 
damped  in  the  vertical  direction  and  extended  in  the 
horizontal  direction.  This  phenomenon  is  being  studied 
in  a  two-dimensional  jet  formed  by  pumping  water  through 
a  rectangular  slit  into  a  channel  filled  with  water.  Hot- 
film  anemometry  techniques  are  being  used  to  measure 
velocity  fluctuations  in  the  jet  and  in  the  immediate 
neighborhood  outside  the  jet.  Since  hot-film  probe 
fouling  cannot  be  avoided  when  making  measurements  in 
water,  a  sensor  calibration  method  which  includes  the 
effect  of  probe  fouling  has  been  devised. 


INTRODUCTION 


When  turbulent  jets  or  wakes  interact  with  an  air-water 
interface,  the  turbulent  eddies  are  apparently  damped  in  the 
vertical  direction  and  extended  in  the  horizontal  direction. 
Experimental  methods  useful  for  studying  turbulent  velocity 
correlations  in  two-dimensional  jets  have  been  devised.  A  unique 
feature  of  this  work  is  the  development  of  a  sensor  calibration 
method  which  reduces  the  labor  required  for  calibrating  hot-film 
probes  used  in  a  multi-channel  array  in  water. 


JET  FLOW  SYSTEM 


A  two-dimensional  jet  is  formed  using  the  flow  system  shown  in 
Figure  1.  Water  is  pumped  upward  from  a  water  channel  0.45  by  0.45 
by  7.5  meters  containing  still  water  using  a  stainless  steel 
centrifugal  pump  driven  by  a  2  HP  motor.  Water  flows  from  the  pump 
discharge  through  polyvinyl  chloride (PVC)  piping,  through  an  inlet 
manifold,  and  into  a  converging  flow  chamber  which  is  submerged  in 
the  channel.  The  water  flows  through  this  converging  section  part 


of  which  is  packed  with  shredded  stainless  steel  and  through  a 
rectangular  slit  formed  by  two  flat  plates  356  mm  wide  by  305  mm 
long  spaced  6mm  apart  and  back  into  the  channel  filled  with  still 
water.  The  flow  rate  through  the  slit  can  be  adjusted  by  adjusting 
the  valves  which  control  the  fraction  of  the  flow  passing  through 
the  by-pass  line.  Since  the  main  circulating  pump  operates  with  a 
suction  lift,  appropriate  valves  and  piping  are  provided  for 
priming  the  pump. 


DATA  ACQUISITION  AND  ANALYSIS  METHODS 


Turbulent  velocity  fluctuations  are  measured  using  hot-film 
probes  connected  to  a  TSI,  Inc.  Model  IFA-100  four-channel 
anemometer  circuit.  By  positioning  four  probes  at  unequal  spacing 
in  the  jet — for  example,  at  0.1,  0.3,  0.7,  and  1.5  meters 
downstream  from  the  slit — data  can  be  acquired  which  will  enable 
calculating  the  spatial  correlations  of  the  turbulent  velocity 
fluctuations  at  six  different  separation  distances — for  example,  at 
separations  of  0.2,  0.4,  0.6,  0.8,  1.2,  and  1.4  meters.  The  analog 
data  are  available  continuously  on  each  of  the  four  channels  of  the 
anemometer  unit.  The  data  acquisition  microcomputer  used  to  record 
the  data — Heathkit  Model  H-8 — samples  the  output  of  the  four 
channel  analog-digital  converter  channel  by  channel  at  2  ms 
intervals, so  there  is  at  least  a  2  ms  lag  between  the  data  from 
different  channels  and  an  8  ms  sampling  interval  for  the  data  on 
any  one  channel.  This  means  that  in  order  to  have  simultaneous 
values  of  the  velocity  fluctuations  to  use  for  calculating  the 
values  of  the  correlation  functions,  some  sort  of  interpolation 
scheme  must  be  used  for  each  channel.  At  any  particular  time 
chosen  for  a  correlation  function  calculation,  there  will  be  one 
actual  data  value  and  three  interpolated  data  values.  The 
calculation  can  be  done  at  2  ms  intervals  using  a  real  data  value 
from  a  different  channel  in  turn  for  each  calculation.  Using  these 
data,  the  wave  number  spectrum  can  be  determined  at  six  different 
wave  numbers  using  a  Fast  Fourier  Transform  method.  If  more  than 
six  wave  number  values  are  desired,  the  correlation  function  must 
be  fitted  to  some  equation  and  interpolated  values  at  intermediate 
separation  distances  must  be  calculated.  An  alternate  method  for 
obtaining  a  spectrum  using  more  wave  numbers  would  be  to  change 
the  probe  spacing  and  collect  more  data.  A  spectrum  calculated 
from  such  data  will  contain  values  calculated  from  velocity 
fluctuation  data  not  all  measured  at  the  same  time. 


PROBE  FOULING  AND  CALIBRATION 


Commercially  available  hot-film  probes  specifically  designed 
for  use  in  water  are  being  used  for  the  measurements.  These  probes 
are  0.152  mm  in  diameter  by  2.0  mm  long  and  are  coated  with  quartz 
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to  minimize  corrosion.  Each  probe  is  heated  to  a  constant 
temperature  when  it  is  immersed  in  the  flowing  water.  The 
temperature  is  set  by  controlling  the  probe  resistance  using  an 
automatic  balancing  Wheatstone  bridge  circuit.  The  sensor  loses 
heat  to  the  surrounding  fluid,  and  the  rate  of  heat  loss  depends  on 
the  heat  transfer  coefficient  for  heat  transfer  from  the  sensor  to 
the  fluid  flowing  past  it.  This  heat  transfer  coefficient  is 
mainly  influenced  by  the  fluid  velocity  perpendicular  to  the 
sensor. 

Since  the  rate  of  heat  loss  is  not  a  linear  function  of  fluid 
velocity,  and  since  there  are  a  number  of  complicating  factors — the 
sensor  supports  influence  the  heat  transfer,  for  example — which 
affect  the  relationship  between  sensor  heat  loss  rate  and  fluid 
velocity,  each  sensor  must  be  calibrated.  The  anemometer  circuit 
output  is  the  voltage  required  to  maintain  the  required  current 
through  the  sensor  in  order  to  maintain  the  sensor  temperature  (  or 
sensor  resistance)  constant.  Sensor  calibration  is  accomplished  by 
immersing  the  probe  in  water  flowing  at  a  known  velocity  and  noting 
the  anemometer  circuit  output.  The  sensor  calibration  system  being 
used  is  simple  to  operate  and  can  be  operated  simultaneously  with 
the  jet  flow  unit  using  the  same  water  that  flows  through  the  slit 
which  forms  the  turbulent  jet.  The  calibration  unit  forms  a 
laminar  jet  in  a  small  enclosed  chamber.  Water  from  the  water 
channel  is  pumped  through  a  calibrated  variable  area  flow  meter  and 
into  the  inlet  of  this  chamber.  The  water  passes  from  the  inlet 
section  through  an  orifice  8.7  mm  in  diameter  into  another  section 
of  the  chamber.  The  water  which  passes  through  the  orifice  forms 
the  laminar  jet,  and  the  velocity  of  this  jet  can  be  varied  by 
varying  the  flow  rate.  The  probe  to  be  calibrated  is  positioned  in 
this  laminar  jet,  and  the  relationship  between  anemometer  output 
and  jet  velocity  is  determined.  The  water  from  the  calibration 
device  passes  out  the  exit  port  and  is  returned  to  the  water 
channel . 

One  annoying  aspect  of  making  turbulence  measurements  in  water 
is  that  the  sensors  become  fouled  quite  rapidly.  This  cannot  be 
avoided,  so  some  method  must  be  devised  for  including  probe  fouling 
in  the  sensor  calibration  scheme.  Since  probe  fouling  always 
occurs  when  hot-film  anemometers  are  operated  in  water,  it  is 
important  that  probe  calibration  be  carried  out  continuously  during 
the  course  of  the  turbulence  measurements.  To  reduce  the  labor  of 
calibrating  several  probes  several  times  during  the  course  of  the 
experiments  carried  out  on  any  particular  day,  is  desireable  to  be 
able  to  calibrate  the  probes  without  moving  or  remounting  them. 

The  objective  of  the  work  was  to  find  out  whether  or  not  the 
probes  foul  in  a  predictable  manner  and  whether  or  not  all  probes 
foul  in  the  same  way  when  operated  in  the  same  water. 

The  results  of  the  calibration  experiments  show  that  probe 
fouling  can  be  expressed  by  the  bridge  output  decreasing  linearly 
with  time.  In  many  cases, the  linear  decrease  is  interrupted  when 
the  water  is  not  moving  relative  to  the  probe.  This  is  illustrated 
in  Figure  2.  At  any  particular  velocity,  each  segment  of  the 


bridge  output  voltage  versus  time  function  can  be  represented  by  a 
linear  function  stored  in  the  mainframe  computer  as  part  of  the 
data  analysis  code.  Such  calibration  curves  are  prepared  during 
the  course  of  turbulence  measurements  with  the  calibration 
extending  over  the  whole  period  during  which  turbulent  velocity 
measurements  are  being  made  in  the  two-dimensional  jet.  The  clock 
time  at  which  any  particular  measurement  is  made  in  the  jet  is 
recorded,  and  the  appropriate  calibration  curve  for  that  particular 
measurement  is  developed  by  reading  the  anemometer  output  voltage 
at  different  velocity  values  at  the  clock  time  noted  for  the 
measurement.  These  data  are  the  fitted  to  a  polynomial  calibration 
equation — 

U  *  Aq  +  AXE  +  A2E2  +  A3E3  +  A4E4 

This  equation  is  then  used  for  interpreting  the  anemometer  output 
data  obtained  with  the  probe  immersed  in  the  turbulent  jet. 

The  calibration  data  obtained  for  several  different  probes 
coincide  closely  enough  so  that  it  is  possible  to  monitor  the 
fouling  of  probes  used  for  turbulence  measurements  in  the  plane  jet 
by  using  a  separate  probe  located  in  the  calibration  chamber.  An 
example  of  the  calibration  data  for  two  different  probes  is  shown 
in  Figure  3.  At  the  start  of  the  calibration  period,  there  is  an 
interval  for  which  the  data  for  the  two  probes  do  not  coincide  very 
closely.  After  this  "seasoning'’  period,  the  calibration  data  do 
coincide. 

The  method  to  be  used  for  incorporating  fouling  into  the  probe 
calibrations  while  avoiding  the  need  to  move  the  probes  during  a 
series  of  turbulence  measurements  is  as  follows. 

1.  Mount  all  probes  in  the  water  channel  and  let 
them  operate  in  the  measuring  mode  for  the 
seasoning  period. 

2.  Install  each  probe  in  turn  in  the  calibration 
chamber  and  obtain  calibration  data  for  each. 

3.  Mount  the  probes  to  be  used  in  the  channel  in 
the  required  positions.  One  probe  to  be  used 
for  monitoring  the  extent  of  fouling  on  all 
probes  is  left  mounted  in  the  calibration 
device.  During  the  course  of  turbulence 
measurements,  calibration  data  are  obtained 
frequently  using  this  monitoring  probe. 

4.  Plot  the  calibration  data  as  E(t)  with  U  as 
a  parameter,  fit  the  calibration  data  with 
linear  functions,  and  use  these  functions  for 
interpreting  the  turbulent  velocity  fluctuation 
data  obtained  from  the  two-dimensional  jet. 

In  order  to  use  this  method  to  include  probe  fouling  in  the 
calibration  functions  used  for  interpreting  turbulence  data,  there 
must  be  at  least  one  actual  calibration  point  for  each  probe.  This 
point  should  be  obtained  after  the  period  of  "seasoning"  in  the 
same  water  as  is  being  used  for  turbulence  experiments. 
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Figure  3.  Comparison  of  Two  Probes 
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The  Skewness  Function  of  Turbulence 


George  Trevino 

Mechanical  Engineering-Engineering  Mechanics  Dept. 
Michigan  Technological  University 
Houghton,  MI  49931 


Abstract 

An  invariant  algebraic  form  of  the  skewness  function  for  a  par¬ 
ticular  type  of  nonhomogeneous  turbulence  is  proposed.  This 
form  is  identical  for  both  the  intended  type  of  turbulence  and 
classical  isotropic  turbulence,  and  as  such  provides  a  unifying 
link  between  turbulence  whose  statistics  are  spatially  varying 


and  turbulence  whose  statistics  are  spatially  constant. 


I.  Introduction 


In  a  recent  publication  (Trevino,  1985),  it  was  suggested  that  traditional 
schemes  for  algebraically  modelling  the  skewness  function  of  turbulence  ty¬ 
pically  lead  to  unfortunate  analytical  difficulties,  especially  when  the 
turbulence  is  nonhomogeneous ;  in  particular  it  was  hinted  that  classical 
methods  do  not  permit  the  absolutely  invariant  structure  of  nonhomogeneous 
turbulence  to  be  explicitly  emphasized  (see  also  Trevino,  1982).  It  was 
therein  also  suggested  that  there  nonetheless  exists  an  invariant  algebraic 
representation  for  the  skewness  function  of  such  turbulence,  but  it  can  be 
formulated  only  for  a  peculiar  class  of  nonhomogeneous  behavior.  That  even 
the  most  general  nonhomogeneous  turbulence  itself  exhibits  a  characteristic 


invariant  structure  somewhere  in  its  hierarchy  of  stochastic  moments  can 

also  be  fundamentally  established  (cf.  Trevino,  1984),  and  how  one  uses  this 

knowledge  is  the  subject  of  a  forthcoming  communication  (Trevino,  to  appear). 

The  purpose  of  this  note  is  to  propose  a  representation  for  the  skewness 

+  is 

function,  S^j^(g,r)  =  <u^(x)Uj(x  +  q)u^(x  +  r)>'  ,  which  identical  for  both 
isotropic  turbulence  and  a  unique  type  of  nonhomogeneous  turbulence.  The 
skewness  function  only  is  to  be  addressed  here,  but  in  principle  any  rele¬ 
vant  three-point  tensor  can  be  similarly  represented.  The  advantage 
provided  by  the  resulting  algebraic  form  of  S_^(q,r)  is  that  it  fur¬ 
nishes  a  unifying  link  between  turbulence  whose  statistics  are  independent 


of  spatial  location  (viz.  x)  and  turbulence  whose  statistics  are  not.  The 


information  available  in  such  a  formulation  is  valuable  in  experimental 


studies  of  turbulence,  where  the  time-decay  of  an  isotropic  nonstationary 
turbulence  is  duly  analyzed  in  terms  of  the  spatial-decay  of  a  nonhomogeneous 


The  C-dependence  of  u  is  here  suppressed  for  convenience. 


stationary  turbulence  (Comte-Bellot  and  Corrsin,  1966);  to  be  sure,  it  re¬ 
quires  that  in  the  advected  frame  of  reference,  the  convective-mechanism  tor 
the  intended  type  of  turbulence  is  characteristically  independent  of  time 
(recall  that  x  -  U  t)  ,  and  that  a  closure  for  same  is  ultimately  obtained 
by  simply  enforcing  the  attendant  invariance  of  • 

2.  The  form  of 

The  skewness  function  is  a  tensor  function  of  the  two  separtion  vectors,  q  and 
r  ;  for  isotropic  turbulence  it  is  written  in  general  as  (cf.  Batchelor,  1967) 

V 

Sin(q,r)  =  I  A  (q  •  q  ,  q  •  r  ,  r  •  r)F^  (q,r)  ,  (1) 

J  n=l  n . '  * 

where  q  •  q  ,  q  •  r  ,  r  •  r  are  the  absolute  scalar  invariants  of  the  full  orthogo¬ 
nal  group .  The  well-known  convective-tensor  of  classical  turbulence  results 
from  Equ.  (1)  by  setting  either  q  =  0  ,  or  q  =  r  ,  to  accordingly  obtain 

sijl(0*£)  *  A(r2)r±r  ^  +  B(r2){ri6jl  +  r  o^}  +  (Kr2)^^  , 

etc..  In  this  note  the  tensor  S, ..(q,r)  is  represented  rather  in  terms  of  the 

ijl  -  - 

invariants,  e^  =  q  •  r/q  •  q  and  -  q  •  r/r  •  r  ,  which  are  unique  to  the  three 
points  specified  by  x  ,  x  +  g  ,  and  x  +  r  ,  and  are  also  more  general  than  the 
standard  metric  invariants  listed  earlier,  e^  and  e^  have  the  advantage  that 
they  remain  invariant  even  though  the  turbulence  itself  may  be  undergoing  some 
scale-change  in  its  statistical  structure,  i.e.  seme  nonhomogeneous  character¬ 
istic  which  can  be  qualitatively  described  by  changes  in  coordinate-axes  scale 
only;  a  nonhomogeneity  of  the  type,  C^(x  +  Ax  ,  r)  =  n  C„(x,£r)  ,  where 
C„(x,r)  =  <u^(x)Uj(x  +  r)>  and  (n,S)  are  simple  scale-factors,  is  a  typical 
example.  In  this  type  of  nonhomogeneity,  (n,?)  are  dependent  only  on  x  ,  Ax 


(and  possibly  t),  but  are  strictly  independent  of  r  .  In  particular,  it 


t  o  I  — 


lows  that  and  e ^  are  such  chat  £q  •  ?r/;q  •  5q  5  ,  etc.,  so  that  these 

invariants  do  not  lose  their  underlying  isotropic  flavor  even  under  the  effects 
of  coordinate-axes  scale-changes .  Note  that  the  classical  isotropic  invari¬ 
ants,  q*q,q*r,r*r,  are  not  preserved  under  scale-change  transformations, 
but  that  e^  and  e^  are  invariant  under  any  combination  of  translations,  rota¬ 
tions,  reflections,  and  the  prescribed  scaling  effects.  This  latter  mode  of 
isotropy  is  more  correctly  denoted  an  affine  isotropy,  in  contradistinction  to 
the  metric  isotropy  of  classical  turbulence  theories. 


The  invariant  form  of  S„^(q,r)  is  herein  to  be  derived  through  successive 
differentiations  of  the  related  first-order  tensor,  T^q.r)  ;  although  not  ger¬ 
mane  to  the  analysis,  this  latter  tensor  could  conceivably  correspond  in  turbu¬ 
lence  theory  to  the  triple-correlation,  <p(x)p(x  +  q)u^(x  +  r)>  ,  where  p(***) 
is  the  fluid  pressure;  accordingly,  the  required  associated  Invariants 
(Batechelor,  loc.  cit.)  of  this  tensor  are 

T^(q,r,a)  =  T  •  a/a  •  a  =  T  •  a 


T(2) (q,r,a)  -  T  •  a/T  •  T  5  T(1) (q,r,a) / j t| 2  ; 


here  a  is  an  arbitrary  unit  vector .  Note  that  both  of  these  forms  are  abso¬ 
lutely  invariant  under  scale-change  transformations  of  the  type  already  con¬ 
sidered,  as  well  as  under  classical  distance-preserving  transformations;  also 
note  that  T^“\***)  is  not  independent  of  T^(***)  .  For  the  special  case 
where  q  =  e  r  both  T^O**)  and  T^(*»»)  are  expressible  in  terms  of  the 
set  of  related  independent  invariants,  { e  ,  r  •  a/r^  ,  r  •  a:  ,  and  the  resulting 
forms  of  T^(e,r)  are 


T^^(e,r)  =  {A(e)  +  r  “B(e)}r\  and  lj”\e,r)  =  {r“A(t)  +  B(e)}  ^r  ; 

A  and  B  are  arbitrary  functions  of  e  only;  because  T^^(***)  and 

are  not  independent,  only  one  of  them,  say  T^(***)  v  T^(***)  ,  need  be 

fully  known  (and  subsequently  used).  The  appropriate  form  of  S„^(q,r)  v 

3T.(e,r)  ,r 

S..,(e,r)  is  now  found  by  determining  —7 — 7 -  ,  and  since  y—  *  r./r  and 

ij  I  ~  dr^jr^  3rj  J 

3e  2 

-  =  -  er./r  it  can  eventuallv  be  established  that 

3rj 

Si^i(e,r)  *  A(e,r)rir^r1  +  BCe.rJr^^  +  C(e,r)r  6^  +  D(e,r)r16i^ 

is  the  intended  form  of  the  skewness  function  of  turbulence.  The  conservation 
of  mass  condition  provides  two  partial  differential  equations  between  the  four 
(4)  functions,  A(***)  ,  B(***)  ,  C(***)  ,  D(***)  ,  so  that  only  two  of  these 
are  independently  arbitrary.  This  same  result  was  originally  reported  by 
Proudman  and  Reid  (1954),  but  the  functions  here  are  not  the  same  as  those  in¬ 
tended  by  them  since  these  apply  to  a  certain  type  of  nonhomogeneous  turbu¬ 
lence  as  well.  The  functions  A,B,...,  can  be  determined  from  experimental 
measurements  by  choosing  r  *  (r,0,0)  (and  e  any  suitable  set  of  numbers)  to 
obtain  the  non-zero  correlations. 


Slll(e,r)  *  Ar  +  (B  +  C  +  D)r  ,  S^0(e,r)  =  s133(e,r)  =  B  r  , 
S,i0(e,r)  *  S313(e,r)  =  C  r  ,  S^U.r)  =  §331(e,r)  =  D  r  . 

3.  Concluding  remarks 

In  dynamic  analyses  the  applicable  forms  of  S„^(e,^)  are  S„^(0,r)  and 
S...(l,r)  ;  these  correspond  to  <u . (x)u . (x) u, (x  +  r)>  and 

*-u.(x)u.(x  +  r)u  (x  +  r )■>  respect ivelv ,  and  Che  number  of  independently  arbi- 

i  -  j  '  -  1  - 


*»’  v"  *..■*  •  *“ . 


trary  functions  needed  to  specify  each  of  these  is  accordingly  one  (1).  These 
forms  are  characteristically  independent  of  x  ,  even  though  the  turbulence  it¬ 
self  may  be  fully  nonhomogeneous .  This  more  general  invariance,  again,  results 
from  the  use  of  related  non-metric  invariants  to  algebraically  represent 
S_i(q,r)  ,  and  is  true  only  for  nonhomogeneous  turbulence  which  undergoes 
coordinate-axes  scale-changes  only.  The  particularly  derived  algebraic  form  of 
is  valid  for  the  collinear  vector  configurations  only  (i.e.  q  =  er)  , 
but  the  method  readily  generalizes  to  non-collinear  vectors  by  introducing  the 

collinear  ratios,  e  =q/r  ,e  =q/r  ,e  =q/r  (Trevino,  1985);  col- 
- x  x  x  y  y  y  z  z  z  - 

linear  ratios,  note,  are  invariant  even  for  scale-changes  which  differ  from 
one  coordinate-axis  to  another,  viz.  £  -*■  (£x,£^,£z)  • 

The  physical  meaning  of  the  above  result  is  that  the  role  of  inertia  forces  in 
the  depicted  type  of  nonhomogeneous  turbulence  (call  it  self-similar  or  self- 
preserving  nonhomogeneous  turbulence),  as  herein  algebraically  expressed,  is 
identical  to  that  for  isotropic  turbulence,  viz.  to  transfer  energy  from  one  part 
of  wave-number  space  to  another  without  changing  the  total  amount  associated 
with  any  particular  directional  component.  In  other  words,  the  related  convec¬ 
tive  terms  do  indeed  affect  the  distribution  of  energy  but  they  do  not  affect 
the  convolution  and/or  modulation  that  due  to  the  nonhomogeneous  nature  ot  the 
turbulence  exists  between  the  energy  associated  with  a  given  wave-number  compo¬ 
nent,  say  d  Z(k)  ,  and  the  energy  associated  with  a  different  wave-number  compo¬ 
nent,  say  d  Z(k')  ;  pressure  terms  and  viscous  terms  are  the  only  terms  that  m 
the  concomitant  hydrodynamic  turbulence  can  influence  this  component  interac¬ 
tion.  Unlike  homogeneous  turbulence,  the  double  correlation  between  different 
wave-number  components  of  nonhomogeneous  turbulence  is  such  that 


<d  Z(<)d  #  0  when  <"  £  < 


i.e.  d  Z(<)  and  d  Z*(<')  are  non-or thogonal  (Trevino,  1981,  1982),  but  the 
triple-correlation,  <dZ(OdZ*(ic')dZ(K~)> ,  is  such  that  dZ(ie)  ,  dZ*(<')  , 
dZ(ic  )  are  orthogonal  for  <  ^  k'  -  ;  *  v  complex  conjugate. 

In  closing  it  remains  exclusively  to  demonstrate  skewness  invariance  for  the 
suggested  turbulence  through  purely  ad^  hoc  geometrical  arguments  on  the  spatial 
evolution  of  the  probability  density  function  (pdf)  itself.  For  the  attendant 
case  of  two  (2)  scaling  effects,  the  most  general  self-preservation  in  the  pdf 
requires  that 


p{u(x  +  Ax)}  =  8  p(au(x)}  , 

where  a  ,  0  are  strictly  ^x-dependent  scale-factors  and  u  is  the  longitudinal 
component  of  u  ,  i.e.  u  =  (u,v,w)  .  In  terms  of  coordinate-axes  transforma¬ 
tions  the  relevant  algebraic  scheme  is 

u(x  +  Ax)  v  u'  =  a  u 

p{u(x  +  Ax)}  v  p'  =  8  p  , 

where  a  3  =  1  ,  i.e.  the  transformation  is  area-preserving ,  this  constraint 
demanded  by  the  pdf  property  that  /  p(5)d£  2  1  always.  Skewness,  recall,  is 
the  f irst  measure  of  dis-symmetry  in  the  pdf  (the  second  measure  is  super 
skewness,  etc.),  and  simple  coordinate-axes  scale-changes  of  the  intended 
type  cannot  alter  this  characteristic  (e.g.  an  odd  function  remains  odd  under 
homogeneous  coordinate-axes  "stretching"  or  "shrinking") .  Skewness  therefore 
is  preserved  in  the  considered  self-similar  evolution.  Corresponding  argu¬ 
ments  apply  to  the  case  where  p  'v  p(v)  ,  p  v  p(v)  ,  p  ^  p(u.v,w)  *  p(u)  , 


etc.,  and  skewness  invariance  is  easily  verified  in  the  special  case  where 
p{u(x)}  is  symmetric  about  the  mean-value  for  all  x  (for  which  the  skewness 
is  always  zero) . 
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Abstract 

The  invariance  concepts  of  affine  geometry  are  invoked  to  estab 
lish  time-independence  for  the  convective  mechanism  of  self¬ 
similar  isotropic  turbulence.  The  decay-law  for  the  integral 
scale  during  the  initial  period  follows  immediately. 


2  s  cf  R  CM,  d 


1.  Introduction 


Affine  geometry  and  affine  transformations  have  recently  begun  to  find  their 
way  into  the  scheme  of  "mechanics-type"  problems  in  such  fields  as  flutter 
analysis1  shaky  structures^1  and  turbulence  their  use  in  qualita¬ 

tively  describing  elastic  deformations  in  solid  mechanics  has  been  known  for 
14 

some  time  .  The  purpose  of  this  note  is  to  adapt  recent  research  re- 
15-19 

suits  ,  regarding  invariant  algebraic  representations  of  nonstationary 
stochastic  processes,  to  the  determination  of  the  characteristic  time- 
invariant  structure  of  self-similar  isotropic  turbulence;  a  closure  for 
same  is  accordingly  formulated,  and  the  decay  of  the  integral  scale  during 
the  initial  period  follows  immediately.  Self-similar  turbulence,  because 
it  can  be  geometrically  described  by  uniform  scale-changes  along  the  ab¬ 
scissa  and  ordinate  axes  respectively,  can  be  characterized  as  turbulence 

20  21 

which  is  invariant  under  the  affine  group  ’ 

The  notion  of  a  self-similar,  or  self-preserving,  turbulence  was  introduced 

22 

by  Karman  and  Howarth  in  1938  .  For  the  most  part  the  self-similar  device 

is  principally  a  mathematical  one,  even  though  the  assumption  itself  is 
quite  valid  in  the  "dissipation  range"  of  the  turbulence  and  in  the  "large- 

scale"  structure  of  the  turbulence;  in  the  dissipation  range  the  self- 

2  2  23 

similar  structure  is  f(r,t)  -  1  -  r  /2X  (during  all  phases  of  the  decay  ), 

while  in  the  large-scale  structure  self-similarity  is  reflected  in  the  time- 

independence  of 

The  essential  idea  of  such  a  flow  is  that  during  the  decay,  the  correlation 
function,  f(r,t)  ,  evolves  in  time,  buts  its  geometric  shape  at  any  time  in¬ 
stant  is  similar  to  its  geometric  shape  at  any  other  time  instant;  i.e.,  the 


integral  scale 


and  the  dissipation  scale, 


Che  time,  c  ,  at  which  Che  moments  are  measured).  In  terms  of  the  geometry 
of  decay  this  invariance  is  a  consequence  of  the  fact  that  any  self-similar 
time-evolution  in  the  probability  density  function  (pdf)  of  the  turbulence 
automatically  preserves  skewness;  i.e.,  any  time-decay  such  that 


p{u(x  ,  t  +  At) }  =  sp{Yu(x,t)}  , 

19 

leaves  the  amount  of  "dis-symmetry"  in  p{.».}  unaltered  .  Note  that  self¬ 
similarity  requires  that  C^r.t)  =  <ui(x,t)u.(x  +  r  ,  t)>  evolve  in  time  in 

2 

a  self-preserving  mode,  or  in  more  practical  terms  that  F(r,t)  =  u  (t)f(r,t) 
undergoes  changes  in  its  coordinate-axes  scales  only,  and  no  changes  in  its 
r-functional  form.  In  particular,  self-similar  decay  requires  that 
F(r  ,  t  +  At)  =  n  F(£r,t)  ,  the  scale-factors  7^  and  ^  being  independent  of 
r  *  (r[  ,  and  dependent  only  upon  t  (and  At)  ;  typically,  n  ^  velocity  scale 
factor  while  C  ^  length  scale  factor. 


15-19 

It  then  follows  that  for  the  considered  turbulence  S^^r.t) 

§  .  (r,a)  ,  and  more  importantly  that  ultimately  (u\t)k(r,t)}  is  also  in- 

h 

dependent  of  t  and  dependent  only  on  r  and  a  ;  here  u  =  cu^u^?  ^  u(t)  , 
and  k(r,t)  is  the  single  arbitrary  function  required  to  completely  specify 


in  the  classical  theory  the  r-dependence  of  S^^(*»»)  .  To  be  sure,  for 


self-similar  turbulence  the  tensor  S^^(r,t)  can  be  written  instead  as 


25 


N 


(n) 


S  (r,a)  =  E  g  (r,a)A.  n(r)  , 
ljl  '  n=i  n  ijl  ' 


where  N  is  here  an  unspecified  finite  number  and  the  g  (•••)  are  invariant 

n 

functions  on  the  affine  (r,t)-space;  the  key  result  in  the  above  is  not  the 

value  of  N  nor  the  unique  form  of  the  g  (•••)  ,  but  rather  it  is  the  explic- 

n 

it  time-independence  of  the  algebraic  representation.  The  equation  for  the 


U 


decay  of  the  integral  scale  is  now  found  by  first  integrating  the  Karrnan- 
Howarth  equation. 


characteristic,  isotropic  turbulence  is  itself  difficult  to  generate  in  prac¬ 
tice,  since  virtually  all  mechanisms  for  creating  turbulence  in  the  laborato¬ 
ry  usually  confer  a  preference  of  at  least  one  particular  spatial  direction 
in  the  flow.  Another  is  that  it  is  common  practice  in  experimental  analyses 
of  turbulence  to  study  the  time-decay  of  an  isotropic  nonstationary  turbu¬ 
lence  by  examining  the  spatial  decay  of  a  nonhomogeneous  (necessarily)  sta¬ 
tionary  (i.e.  steady-state)  turbulence,  and  accordingly  imposing  the  "frozen 
turbulence"  assumption  to  transform  time-derivatives  into  spatial  derivatives 
(and  vice-versa)  as^ 

{— }c  =  U  —  ,  U  'v  mean-flow  velocity. 

3t  following  the  3x  7 

mean  motion 

This  scheme,  although  appealing  and  quite  practical,  has  the  defect  that  in 
a  nonhomogeneous  turbulence  pressure  terms  play  an  important  role  in  the  de¬ 
cay,  even  for  the  total-energy  case;  indeed,  if  the  spatial  evolution  of  the 

covariance  of  such  a  flow  is  in  fact  of  a  self-similar  nature,  then  theoret- 

12 

ical  results  spatially  analogous  to  the  foregoing  establish  that  convec¬ 
tive  terms  are  not  at  all  the  proponents  of  the  "spatial  decay"  of  the  tur¬ 
bulence - pressure  terms  are  the  only  carriers  of  nonhomogeneous  energy- 

transfer  while  convective-transfer  is  necessarily  spatially  constant,  i.e. 
in  a  truly  nonhomogeneous  self-similar  flow  pressure  terms  cannot  be  ne¬ 
glected.  In  the  study  of  total-energy  decay  of  homogeneous  isotropic  tur- 

23 

bulence  pressure  terms  are  not  only  negligible,  they  are  indentically  zero 

Decay  laws  of  the  above  form  have  also  not  otherwise  been  reported  in  the 

27 

literature  (see,  for  example,  P.G.  Saffman  );  in  fact,  in  order  to  obtain 
with  the  herein  proposed  analysis  decay  results  comparable  to  Saffman’ s,  n 
in  the  derived  power-law  decay  for  kinetic  energy  would  have  to  be  n  *  -3/5  . 


6 


Saffman's  results  are  obtained  by  imposing  the  analytical  consequences  of 
the  invariance  of  a  particular  integral  of  the  velocity-covariance  tensor, 
C^(r,t)  and  not  by  imposing,  as  done  here,  the  consequences  of  the  newly 
found  time-invariance  of  the  inertial-transfer  tensor. 


In  closing  it  remains  to  iterate  that  the  above  results  can  only  be  estab¬ 
lished  through  affine  invariants  modelling,  i.e.  by  modelling  using  parame¬ 
ters  that  remain  invariant  even  though  the  turbulence  itself  is  undergoing 

temporal  changes  in  its  coordinate-axes  scales - it  cannot  be  established 

through  metric  invariants  modelling  without  introducing  some  germane  postu¬ 
late  about  the  metric  structure  of  .  In  other  words,  by  modelling 

the  related  self-similar  structure  of  S^^(***)  through  non-metric  invariants, 
invariants  which  are  unique  and  peculiar  to  this  type  of  evolution,  it  is 


then  possible  to  analytically  establish  a  rather  fundamental  characteristic 


about  the  convective-transfer  in  self-similar  turbulent  decay. 
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